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1. Introduction 

As numerical simulations in lattice QCD mature, we can attack more difficult or complicated 
problems in strong interactions. Calculations of hadronic matrix elements of electroweak operators 
are predictions from lattice QCD, which can be used to determine parameters of the standard model. 
Another interesting application is to investigate interactions among hadrons using lattice QCD. 
Since hadrons are bound states of quarks and gluons, their interactions are residual effects of strong 
interactions among them inside hadrons, and therefore they are more complicated quantities in 
lattice QCD than properties of an isolated hadron such as masses or decay constants. In addition, 
the euclidean nature of lattice QCD makes problems of hadron interactions more complicated and 
difficult. Thanks to the finite volume technique [jjj], however, scattering lengths and phase shifts of 
hadrons can in principle be extracted from lattice QCD. Some recent progress on these quantities 
is reviewed in this paper. 

I also consider a more complicated but important quantity of hadron interactions, the force 
between nucleons (the nuclear force). In 1935 Yukawa introduced a virtual particle, the pion, to 
account for the nuclear force[Q], by which protons and neutrons are confined in nuclei. Since 
then enormous efforts have been made to understand the nucleon-nucleon (AW) interaction at low 
energies both theoretically and experimentally. In Fig. Ill, I present modern AW potentials, which 
are characterized by the following features^, |J]. At long distances (r > 2 fm ) there exists weak 
attraction, which is well understood and is dominated by the one pion exchange (OPE), as first 
pointed out by Yukawa. At medium distances ( 1 fm < r < 2 fm), contributions from the exchange 
of multi-pions and heavy mesons such as p, CO and a lead to slightly stronger attraction. At short 
distances (r < 1 fm), attraction turns into repulsion, and it becomes stronger and stronger as r 
gets smaller, forming the strong repulsive core[||], which is essential not only for describing the 
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Figure 1: Three examples of the modern AW potential in 'So (spin singlet and 5-wave) channel: Bonn[|^], 
Reid93[§ andAV18[0. 
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Figure 2: 5-wave scattering lengths in unit of 1 jM n . 



NN scattering data, but also for the stability and saturation of atomic nuclei, for determining the 
maximum mass of neutron stars, and for igniting Type II supernova explosions [||]. Although the 
origin of the repulsive core must be related to the quark-gluon structure of nucleons, it remains one 
of the most fundamental problems in nuclear physics for a long time[j7|]. It is a great challenge for 
us to derive the nuclear potential including the repulsive core from lattice QCD. Recent progress 
on this issue is explained in this paper. 



2. Conventional method in lattice QCD 

In this section, some recent results for the scattering length and phase shift for pions are dis- 
cussed. 



2.1 Scattering length of pions 

In Fig. [2], recent theoretical estimates of the S-wave %% scattering lengths, and c?q in units 



of 1 /M K , are summarized [11], where L and / of the scattering length a£ represent the total angular 
momentum and the total isospin, respectively. The band bounded by two black solid lines repre- 
sents the theoretically allowed region (universal band) [|lT|] . The three black dots are predictions 
from 2-flavor chiral perturbation theory (ChPT) at tree, 1-loop and 2-loop levels from left to right. 



Other results are obtained by using the formulae! 12] 
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where M^C^ are corrections to ChPT from dispersion integrals. They are approximated at the 
physical point as 

M*a 3 = 0.135 + 0.77(a§- 0.220) -1.50(a§ + 0.0444), (2.3) 
M*a 4 = 0.061 +0.48(ag- 0.220) -0.26(4 + 0.0444). (2.4) 

^4 are the low energy constants in Nf = 2 ChPT and defined by 

M 2 = M 2 {l + ^l 3 + 0(x 2 )}, F K =F{l+xk + 0{x 2 )} (2.5) 

with M 2 = 2Bm and x = M 2 / ( 16tt 2 F 2 ). Here m is the quark mass and F is the pion decay constant 
in the chiral limit. Once % and I4 are known, the above expressions give Oq and a^. 

The red ellipse in the figure corresponds to I3 = 2.9(2.4) from the mass spectrum of the pseu- 
doscalar octet and l 4 = 4.4(2) from the scalar form factor of the pion. The narrow azure strip 
indicates the region allowed if T3 is treated as free parameter while I4 is fixed to the above value. 
There exist several lattice determinations of It, and I4. The MILC collaboration[]l3|], using 2+1 fla- 
vors of staggered dynamical quarks, obtain the low energy constants £4 5,6,8 of Nf = 3 ChPT, which 
are translated to I3 = 0.8(2.3) and I4 = 4.0(6) by standard one loop formulae. The results for a[J 
and a 2 , are given by a green ellipse in the figure, which is consistent with the previous estimate (the 
red ellipse). The result of Nf = 2 dynamical Wilson quarks, T3 = 3.0(5) (l)[jl4"|], with free I4 leads 
to the narrow purple strip in the figure, which intersects the red ellipse. 



Recent results of the ETM collaboration[|5]], h = 3.65(12) and l 4 = 4.52(6), give M n a 



0.221(1) and M n a^ = -0.0427(1) (with M K = 140 MeV and F % = 93 MeV), which nicely agree 
with the red ellipse with smaller errors. 

These lattice results for a{] and a 2 , are indirectly obtained by using Ij and l 4 , which have been 
determined from the quark mass dependences of M n and F K . There exist several direct estimates 
on <3q using the finite volume method[[l]]. The horizontal orange band in Fig. || indicates the direct 



result for by the NPLQCD collaboration [[16J] with domain-wall valence quarks on a 2+1 stag- 
gered sea. Recently they have reported a more accurate result, M n a^ = —0.0443 (4) Jl7|] , which is 
consistent with the red ellipse. So far ChPT predictions and all indirect /direct lattice results agree 
with each other. However a direct evaluation of a[J is missing. Thus it is a great challenge for lattice 
QCD to extract a[J directly by the finite volume method. 

2.2 The p meson decay width 

A less difficult than £?{j but still challenging problem is the determination of the P-wave scat- 
tering phase shift for the 1=1 two-pion system, from which the p meson decay width is extracted. 
Applying the finite volume method in the laboratory system to this problem, the CP-PACS collab- 
oration has calculated the decay width[jl8|] in lattice QCD with a renormalization group improved 
gauge action and Nf = 2 dynamical clover quarks on a 12 3 x 24 lattice at m K /m p ~ 0.41 and the lat- 
tice spacing 1/a ~ 0.92 GeV. In order to realize the kinematics such that the energy of the two pions 
is close to m p , one pion has a non-zero momentum p = (2n/L)e.T, and the other is at rest, while the 
p meson has the same non-zero momentum. Energies of these states for non-interacting hadrons 
are W® = y/ m 2 + p 2 + m % for the two pions and W 2 ° = a/ m p + P 2 f° r tne P- At m^/ntp ~ 0.41 on 
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Figure 3: Scattering phase shift sin 2 5(lower panel), positions of m p and resonance mass Mr (upper panel) 
in lattice units. "Cont" refers to results obtained with the continuum dispersion relation while "Lat" to those 
obtained with the lattice dispersion relation. 



a 12 3 x 24 lattice, the invariant mass of two free pions takes the value y/s ~ 0.97 x m p , which is 
much closer to m p than that in the center of mass system, E = 2\/m 2 + p 2 ~ l.47m p . The hadron 
interaction shifts the energies from W® to W n , which are related to the two-pion scattering phase 



shift 8 in the infinite volume through the Rummukainen-Gottlieb formula[19]. 

To extract two energy levels W n close to each other, a 2 x 2 matrix of the time correlation 
function 

= I <o|(*w)t(0(*w)(f,)|o> (0|(jw)t(0fl3(O|0> \ r2 fi 

1 ' [ (o\pl(t)(7in)(t s )\o) <o|/>J(Op3(0|o> ) K ' ) 

has been constructed, where p?,{t) is an interpolating operator for the neutral p meson with the 
momentum p and the polarization vector parallel to p, and (im)(t) is an interpolating operator for 
two pions, (%%)(t) = (7T~(p,?)7r + (0,?) — 7r + (p,f)7r _ (0,f))/v / 2- Two energy levels are obtained 
from two eigenvalues X„(t,t R ) of the matrix M(t,t R ) = G{t)G~ l {t R ) as 

Ut,tR)=e- W " {t - tR) (2.7) 

for large t. To keep symmetry between source and sink in G(t) U(l) noises are introduced. The 
total number of quark propagators per configuration becomes 520, including 10 noises times 2 



source points J 18|] 



From the energy levels W n extracted, the invariant mass and the momentum k are given by 
yfi = y/W^-p 2 and k 2 = s/4 — m 2 , assuming the continuum dispersion relation, while they be- 
come cosh(- v /s) = cosh(W„) — 2sin 2 (/?/2) and 2sin 2 (&/2) = cosh(y / s/2) — cosh(m^), using a lat- 
tice dispersion relation. Then W\ leads to tan 8 = 0.0791(9) at = 0.776(8) (continuum disper- 
sion) or tan<5 = 0.074(7) at y/s = 0.800(8) (lattice dispersion), and W 2 gives tan<5 = -0.38(11) 
at y/s = 0.95(3) (continuum) or tan 8 = -0.52(14) at y/s = 0.98(3) (lattice). In both cases it is 
observed that 8 > (attractive) at y/s <m p = 0.86(1) while 8 < (repulsive) at y/s > m p . This 
property confirms the existence of a resonance at a mass around m p [|2Q|]. 
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In Fig. [|, sin 2 8, which is proportional to the scattering cross section of the two-pion system, 
is plotted against the invariant mass i/s. The two lines in the figure represent fits of the data by 

tang = g p e 

6n y/s(M%-s) 

where M R is the resonance mass, plotted also in the upper panel of Fig. |3[ and g pK7[ is the effective 
coupling of p%%. If gpmt were obtained at several quark masses, it could be extrapolated to the 
physical quark mass. Since this calculation is performed only at one quark mass, however, it 
is assumed that g p7[K is independent of the quark mass, so that the p meson decay width at the 
physical quark mass is estimated by 

^^S^*'* 4 - 128 ^ <2 - 9) 

where (F h ) 2 = (mf) 2 /4-(mf) 2 . The fit result of g pKS gives TP h = 162(35) MeV (continuum) 
and r ph = 140(27)MeV (lattice), which are consistent with the experiment value, T = 150 MeV. 

Since the first attempt is successful, the next step is the chiral extrapolation. The p meson 
decay width is a good bench mark quantity which shows the dynamical nature of QCD vividly. 
Therefore collaboration groups which have full QCD configurations at small quark masses are 
encouraged to calculate this quantity. 



3. Potentials for heavy hadrons 

In the following two sections, extractions of the potential between hadrons from lattice QCD 
are discussed. A straight-forward way to define a potential is to calculate an energy of a two- 
hadron system as a function of the distance between the two hadrons. A difficulty exists in the 
definition of the distance between two hadrons, since the two hadrons are moving around changing 
their relative distance. To overcome this difficulty, one (infinitely) heavy quark may be introduced 
in each hadron, so that the distance between two hadrons is defined by the distance between the 
two heavy static quarks, which do not move in space for all the time. This definition, similar to 
the static quark potential from a Wilson loop or Polyakov lines, is indeed employed to calculate 
potentials between two heavy hadrons in lattice QCD. 
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Figure 4: The binding energy, Vbb — 2Vb (GeV), as a function of r (fm). Left: All the light-quark flavors 
are different. Middle: One pair of quark flavors is identical. Right: Two pairs are identical. 
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Figure 5: The central meson-meson potential (MeV) as a function of r in lattice units. 
3.1 Baryon-Baryon 



In Ref. pip, the energy of two heavy baryons has been measured as a function of the distance 
r between two heavy quarks Q, using quenched QCD with the plaquette gauge action and the 
Wilson quark action at a ~ 0.19 fm on a 20 3 x 24 lattice. The binding energy of two baryons, 
Vbb — 2Vb, is plotted as a function of r in Fig. [l] at m n ~ 500 MeV for the light quark, where Vbb is 
an energy for two baryons while Vb is an energy for one baryon. The differences of the three figures 
are explained in the caption. Surprisingly, the binding energy is very small and shows almost no 
dependence on r for all three cases. Thus, a repulsive core is not seen for the heavy baryon potential 
in this calculation. 

3.2 Meson-Meson and others 



In Ref. p2p, the binding energy of two heavy-light mesons has been computed, each of which 
is made of one static quark and one light quark, using quenched QCD with the DBW2 gauge action 
at a ~ 0.1 fm and the Wilson quark action at m n ~ 400 MeV. The central potential for total spin S 
and isospin / of two light quarks in a meson-meson system is decomposed as 



Vi,s{r) = Vi (r) + (a! • a 2 )V a {r) + fa • T 2 )V T (r) + {a x ■ C7 2 )(t 1 • T 2 )V ax (r) 



(3.1) 



where (7; (T,) acts on the spin (isospin) of the light quark qi, and r is the distance between the two 
static quarks. Results for Vx (r) (X = 1, a, X, ox) are plotted in Fig. |5[ All potentials show attraction 
at short distance (r < 0.2 fm ). 



In Ref. [23], a meson-meson potential and a meson-baryon potential has been calculated with 
the same setup of the baryon-baryon case in the previous subsection. As shown in Fig. ^, however, 
the dependence on r is very small in both cases. In particular, no short-distance attraction is ob- 
served for the two heavy-light mesons, contrary to the results in Fig. |5[ There are several candidates 
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Figure 6: The meson-meson potential (left) and the meson-baryon potential as a function of r at m-n ~ 700 
MeV. Triangles (circles) denote results without (with) light quark exchange diagrams. 



which might cause the difference: the light quark mass, the lattice spacing, the lattice volume, the 
statistics or the way to extract the binding energy. Further investigations are necessary for definite 
conclusions about the static quark approach to potentials between hadrons. 

4. Potentials from wave functions 



In Ref.[|24|], the nucleon-nucleon (AW) potential has been extracted in lattice QCD by a totally 



new and different method, which is explained in this section. 

4.1 Wave functions and potentials 

The starting point of the new method for the AW potential is a wave function defined by 

<p(r,E) = (0\N(x,0)N(y,0)\W;E), r = x-y, (4.1) 

where N(x,t) is an interpolating field of the nucleon, \2N;E) is a 2N state with energy E < inelastic 
with the inelastic threshold energy inelastic of the 2N system. This wave function was first analyzed 



for a spin model in Ref. Q25D, and has been employed in lattice QCD to calculate the %% scattering 



length[26] 



This wave function has the following properties. For large enough r = |r|, (Ho — E)(p(r,E) = 
-V 2 

0, where Hq = — — is a free non-relativistic Hamiltonian with the reduced mass /J. = mjv/2. This 

means that the interaction between two nucleons vanishes for a large separation. Moreover, for 
example in the case of the 1 Sq channel, it can be shown that 

^(r,£)^^W Sin( ^ +W) +--, (4.2) 
kr 

where So(k) is the S wave scattering phase shift and k is determined by E = k 2 /(2jx). In Ref. [26, 



27], it is shown for a two-pion system that 

ik tjii i\ ■ n \ , n f di P H (P' k )Mpr) 



«p 5 (r,£) = e^ + ^H(t,k)j (kr)+P f * P H{P f J [f } + W s (r), (4.3) 
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where H(p,k) = (E k +E p )M(p,k)/(SE k E p ) with E p = yjm\ + p 2 and M(p,k) is the (off-shell) 
scattering amplitude. A similar but a little more complicated expression can be shown also for the 
AW system[|2j|]. For large r, if E < inelastic > the contribution from /others vanishes exponentially 



and eq.(|4.3|) approaches eq.(|4.2|). 



The potential is defined from the wave function through the Schrodinger equation 

[H + V(r)}(p(r,E) = Eq>(r,E), (4.4) 

which symbolically gives 

v{r) _ 

4.2 Lattice calculation 

A corresponding lattice wave function in a finite box at E ~ is given explicitly by 

<P( r > £ ) = I j\L^ a p(mU^]+^0)N(x,0)\2N;E), (4.6) 



X 



where N' a = £ a bc(' q a CY5T 2 q b )q l a , P 1 = 1 and P° = T 2 are isospin projections, P [ = 1 and P° = 02 
are spin projections, and the summation over x gives zero total momentum. The summation over 
the discrete rotations ffl of the cubic group implies that the state belongs to the A\ representation 
of the cubic group, which is expected to couple to an L = ground state as well as L > 4 excited 
states of the rotation group in the continuum theory. 

The wave function without projections has been extracted from the 4-point nucleon correlator 

as 

F NN (x,y,f;? ) = (0\NUx,t)N J p (y,t) J? NN \0) = £A n (0|<(x,f)^ (y,?)|2A r ;^)^ £ " (f -' o) (4.7) 

n 

where A„ = (2A r ;£ , „| ( ^ NN (f )|0) and c /nn(^o) = P-jP^ptA^vVj . In Jf , a wall source is employed 
by replacing <?(x,?o) with Q(to) = Lx<?( x ) Z o) a f ter Coulomb gauge fixing. For large enough t the 
wave function for E = Eq is obtained. 

In a finite volume, energy levels are shifted from those in the infinite volume as AE n (L) = 
E n (L) —E n (L = 00) = 0(L~ 3 ) due to interactions. From these shifts, the scattering phase shift can 
be extracted pf]. 

4.3 Numerical simulations 

In the actual calculation, the NN scattering for L = is considered. There are two channels, 
the spin singlet (S = 0) channel So and the spint triplet (5=1) channel Si, where the standard 
notation 2S+l Lj is used. Since only a central potential appears for l So, the definition ( |4.5| ) directly 
gives Vnn(i") = Vc(r). On the other hand, in the case of 3 S\, the potential becomes Vnn(i") = 
Vc(r) + Vr(r)5i2(r), where Si2(r) = 3(<Ti • r)(<72 ■ r ) _ (ci • 02) is the tensor operator. The tensor 
potential induces mixing between 3 Si and 3 D\ states, so that the definition (1.5) gives the so-called 
effective central potential, which includes the mixing effect from the 3 D\ state as a second order 
perturbation. 
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Figure 7: Figures from [[24]]. Left: The normalized AW wave function for the singlet(circles) and the 
triplet( triangles) as a function of r (fm) at m K ~ 527 MeV. Right: The central (effective central) AW potential 
(MeV) as a function of r (fm) for the singlet (triplet) at m % ~ 527 MeV. 

Numerical simulations have been made in quenched QCD on a 32 4 lattice with the plaquette 
gauge action and the Wilson quark action, at a ~ 0.137 fm from the p meson mass. The number of 
configurations is 2000 for m n ~ 370 MeV and 527 MeV, and 1000 for «% ~ 732 MeV. A Dirichlet 
boundary condition (DBC) in time and a periodic B.C. in space are employed, and the wall source 
with Coulomb gauge fixing is placed at to = 5, to avoid an influence of the DBC. Calculations have 
been performed on a Blue Gene/L at KEK, which has 57.3 TFlops peak performance. It took about 
4000 hours of 512 Nodes (a half -rack, 2.87TFlops peak) with 34-48% sustained speed to complete 
them. 

4.4 Results 

In the left panel of Fig. [7[ the AW wave function normalized to 1 at r ~ 2.2 fm is plotted 
as a function of r for 1 Sq and 3 Si at m n ~ 527 MeV. In both cases the wave function shows an 
increase at 0.5 fm < r < 1.5 fm, suggesting attraction, while it decreases at r < 0.5 fm, indicating 
the existence of the repulsive core. In the right panel of Fig. [7] the central (effective central) AW 
potential extracted from the wave function by eq.(^^) at m n ~ 527 MeV is plotted as a function 
of r for So ( 3 S\). Interestingly the potential obtained qualitatively agrees with the AW potential 
determined from scattering experiments: weak attraction at long distance, a little stronger attraction 
at intermediate distance and strong repulsion at short distance (the repulsive core). The solid line 
is the Yukawa potential (One Pion Exchange Potential) given by 

^Yukawa / ) = §In [ g) ( g l [ <g) / Y (4 §) 

1 4k 3 \2niN J r ' 

which agrees well with the data at long distance, with g^/(^Jt) c± 14.0 from experiments, m n ~ 
0.527 GeV and ~ 1.34 GeV from lattice data. Thus, in some sense, the Yukawa theory for the 
nuclear force at r > 1 fm is confirmed by lattice QCD. 

The quark mass dependence of the AW potential for 'So is shown in Fig. [| As the quark 
mass decreases, the repulsive core at short distance gets stronger, and at the same time, attraction 



at intermediate distances also becomes a little stronger [29]. In the future it will be interesting to 
see if this quark mass dependence remains in full QCD. 
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5. Theoretical considerations and future perspectives 

5.1 Energy dependence of the potential 

Since the wave function depends on the energy E = p 2 / (2ju), the potential defined by eq.(1.5) 
is also energy dependent in general: 



V J (r,E) 



(E-H )(p J (r,E) 
(p J (r,E) 



(5.1) 



where /, the total angular momentum, is fixed, and r = |r|. To make our argument simpler, the 
spin degrees of freedom are not considered here. Since V J (r,E) carries more information than the 
scattering phase shift 8(E) does, V J (r,E) is redundant and therefore not physical. For the potential 
defined from the wave function to be physically meaningful, its energy dependence must be weak 
for some range of small E. 

In Fig. ^, the %% wave function and the corresponding potential in the 5 channel are shown 
at r = (x,y,0) in quenched QCD at m n /m p ~ 0.51 for p = (left) and p = (2n/L,0,0) (right) on 
an L 3 box[f30"|] . Although the wave functions at the two different energies are very different, the 
potentials look similar. For the Tin system, the energy dependence of the potential is not so strong at 
low energy. Note, however, that the wave function (and therefore the potential) at p = (2jt/L,0,0) 
is calculated at equal time in the laboratory system and is transformed back to the center of mass 
system, so that the relative time between the two pions is x dependent: t = vyx. 

5.2 A unique local potential 



In this subsection a method to extract a unique potential from wave functions is proposed! 28]. 
As mentioned in the previous subsection, the potential seems E dependent in general. In addition, it 
may also depend on the choice of the interpolating operator N(x,t) in (fO|): one can use a different 
N(x,t) unless it changes the asymptotic behavior (the phase shift). 
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Figure 9: Upper plots: The nn wave function at p x = 0,f = (left) and p x = 2n/L,t = vyx (right) in 
quenched QCD at m % jmp ~ 0.51. Lower plots: The corresponding potential. 

To begin with, it is argued that, by introducing a non-local potential, the energy dependence 
of the potential can be removed. The non-local potential is defined by the following equation, 

K(r,E) = (E-H )(p(r,E) = J £u(r/Mr',E), (5.2) 

r> 

where all quantities are taken to be real. For notational simplicity, £ is used even for the case that 
a variable is continuous and the superscript / for the angular momentum is omitted here. Since 
{(p(r,E)}E form a complete set, there exists an inverse such that 

J^q>(r,E)q>- l (Ey) = S(r-r'), ^\E,r) = ^"^i?') (5.3) 

E E' 

where r\ E l E , is the inverse of v\e,e> = Hr 9( r ^E)^)(r,E'). Using this we obtain 

U(r,r')= J £K(r,E)<p- l (E,r'). (5.4) 

E 

In the actual simulations, it is impossible to obtain (p(r,E) for all E. If <p(r,£')'s are obtained 
for E = Eo,E\ ,E n , we expand U in terms of a derivative d such that 

U{rJ) = [U (r)+U l (r)d + U 2 (r)d 2 + ---]8(r-r') = 1 £U k (r)d k 8(r-r'). (5.5) 

k=o 

The coefficient Ut can be obtained as 

U k (r) = £^(r,£y)^(r) (5.6) 



12 



Hadmn interactions from lattice QCD 



where is the inverse of (j>j^ = d k ([>(r,Ej). Note that d becomes a difference on a lattice. 

It is now argued that the non-local potential U obtained above can be transformed to a local E 
independent potential V. The Schrodinger equations for them become 

(H + U)\(p,E) =E\<p,E), (r\U\/) = U(rS), (r\(p,E) = <p(r,E) (5.7) 
(H Q + V)\ ¥ ,E) = E\y,E), (r\V\/) = V(r)8(r- /), (r\y,E) = y^E), (5.8) 

where Hq + U is not hermitian unless Ufar*) = £/(/,?•), while Ho + V is hermitian. Note here 
that the eigenvalues E in both equations must be equal to ensure that both wave functions give the 
same scattering phase shift (and the bound state spectrum if any). Existence and uniqueness of V 
are suggested by the inverse scattering theory, which tells us that a local potential for a fixed J is 
uniquely reconstructed from the scattering phase shift for all E and informations of possible bound 
states. Since the eigenfunctions \y,E) of the hermitian operator Hq + V satisfy (y,E\y,E') = 
C\8e,e', in contrast to ((p,E\(p,E') = rj EE >, the transformation function (r|A|r ; ) = A(r,r') defined 
by 



y,E) = A\<p,E), \<p,E) = A" 1 1 !//•,£}, (5.9) 



can be constructed as 



A = £ |v,£)T7^,(<p,£'|, A' 1 =Y,\Q,E)-^(Y,E\. (5.10) 

E,E' E L £ 

The compatibility of non-local and local potentials leads to 

A{H Q + U)\<p,E)=EA\<p,E) = (H Q + V)A\(p,E), (5.11) 
which, together with the completeness of \(p,E), implies 

[A,H ] +AU = VA. (5.12) 
If the number of degrees of freedom for E and r is denoted as #E = #r = N on a finite lattice, 



eq.(5.12) gives N x N constraints, while the unkown A and V have N x N and ,/V components, 



respectively. Therefore the determination of both A and V from eq.(5.12) seems an ill-posed prob 



lem. However, N(=#E) components of A can be taken freely, since the eigenvalue equation (| 



does not depend on Ce, the norm of the eigenfunction, for all £"s. Using this freedom, eq. Q5.12D is 



enough to determine A and V. Eq.( |5.12| ) is the master equation, which determines the unique local 
energy independent potential V, once a non-local potential U is obtained from a wave function 
defined through a particular choice of operators. This conceptually solves the uniqueness (or E 
dependence) problem of "the potential from a wave function". The results in the previous section 
correspond to the 0th order solution of this equation: V = Uq and A = 1. 

5.3 Future perspectives 

To understand interaction properties of hyperons, baryons which include at least one strange 
quark, is an important subjects in the nuclear physics. Hyperon-nucleon (YN) and hyperon-hyperon 
(YY) interactions are relevant to structures of the neutron-star core and the existence/absence of H- 
dibaryon states. Properties of hypernuclei, nuclei which contain hyperons, will be also studied, as 
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Figure 10: Figures from pT|]. Left: The normalized wave function of pZ° in the l So (circles) and 3 S\ 
(triangles) channels as a function of r (fm). Right: The corresponding effective central potential (MeV) . 



a project at J-PARC (Japan Proton Accelerator Research Complex). However, YN and YY interac- 
tions are poorly known both theoretically and experimentally so far. Lattice QCD calculations of 
the AW potential in the previous section can be extended to YN and YY potentials. 

In Fig. [T(| the wave function and the corresponding potential for pE are plotted as a function 
of r in quenched QCD [31 1. The lattice parameters are the same as in the case of the AW potential 
in the previous section. The light quark mass corresponds to m % ~ 370 MeV while the strange 
quark mass is tuned to reproduce the physical K meson mass, — 550 MeV. Qualitative features 
of the YN potential are similar to those of the AW potential. Weak attraction appears at long and 
intermediate distances while a strong repulsive core shows up at short distance. However the spin 
dependence of the potential is stronger than in the AW case. In particular, the repulsive core in 
the So channel is much stronger than that in the 3 Si channel. Although these results are still 
preliminary, they are interesting and encouraging as a first step. 

Lattice QCD calculations of potentials between hadrons from wave functions have just begun 
and the first result of the AW potential surprisingly reproduces all the known features of the AW 
potential such as weak attraction at long distance, a little stronger attraction at intermediate distance 
and a strong repulsive core at short distance. For a quantitative comparison between lattice results 
and experimental ones, however, chiral and continuum extrapolations are necessary to remove 
systematic errors. The inclusion of dynamical quark effects in the AW potential will be the most 
exciting improvement in future calculations. Now a door is open for us to vast fields in nuclear 
physics with lattice QCD. 
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